ELSEVIER 


Available  online  at  www.sciencedirect.com 


Journal  of  Power  Sources  152  (2005)  60-66 

www.elsevier.com /locate /jpowsour 


Simulation  of  the  current  density  distribution  for  a  PEMFC  by  using 
measured  electrochemical  and  physical  properties  of  the  membrane^ 

Takuto  Araki* *,  Hironori  Koori,  Takuya  Taniuchi,  Kazuo  Onda 

Department  of  Electrical  and  Electronic  Engineering,  Toyohashi  University  of  Technology  1-1  Hibarigaoka, 

Tenpaku,  Toyohashi,  Aichi  441-8580,  Japan 

Received  30  November  2004;  accepted  27  December  2004 
Available  online  1  August  2005 


Abstract 

In  order  to  grasp  properly  PEFC  power  generation  performances,  it  is  necessary  to  know  aspects  of  water  management  such  as  transmissivity 
and  electro-osmotic  coefficient  of  water  vapor  through  the  membrane,  and  factors  for  power  loss  such  as  active  and  resistive  overpotentials. 
In  this  study  we  have  measured  these  factors  to  analyze  our  experimental  results  of  PEFC  power  generation  tests  by  our  two-dimensional 
simulation  code.  It  considers  simultaneously  the  mass,  charge  and  energy  conservation  equations,  and  the  equivalent  electric-circuit  for  PEFC 
to  give  numerical  distributions  of  hydrogen/oxygen  concentrations,  current  density,  and  gas/cell-component  temperatures.  The  numerical 
distributions  of  current  density  under  various  operating  conditions  agreed  well  with  the  measured  distributions  by  segmented  electrodes, 
which  had  grooves  for  hydrogen/oxygen  supply  and  were  molded  in  our  test  cell  being  electrically  insulated.  Hydrogen/oxygen  concentration 
changes  measured  by  gas  chromatography  along  the  gas  supply  grooves  gave  also  the  experimental  current  distributions,  which  coincided 
almost  with  those  by  the  segmented  electrodes.  Factors  to  correct  the  small  difference  between  the  measured  and  the  calculated  are  also 
discussed  from  the  stand  point  of  the  physical  meaning  of  the  calculated  results  considering  factors  which  are  not  taken  into  account  in  our 
code. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Establishing  the  precise  numerical  model  for  the  proton 
exchange  membrane  fuel  cell  (PEMFC)  provides  a  useful 
method  to  estimate  and  improve  PEMFC  performance.  The 
first  objective  of  the  present  investigation  was  to  build  a 
numerical  model,  which  can  describe  the  PEMFC  perfor¬ 
mance  such  as  a  V-i  characteristics  and  current  distributions 
inside  the  cell.  For  properly  building  the  PEMFC  model,  it 
is  necessary  to  know  the  water  management  factors  such 
as  transmissivity  and  electro-osmotic  coefficient  of  water 
vapor  through  the  membrane-electrode  assembly  (MEA),  and 
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power  loss  factors  such  as  activation  and  resistive  overpo¬ 
tentials,  since  these  factors  have  significant  impact  on  the 
performance  of  the  PEMFC.  Nevertheless,  the  effects  of  these 
factors  on  the  PEMFC  have  not  been  defined  sufficiently. 
Then,  obtaining  the  basic  information  for  PEMFC  is  a  sec¬ 
ond  objective  of  this  study.  A  comparison  of  measured  and 
calculated  current  distributions  is  finally  conducted  to  ver¬ 
ify  the  reliability  of  our  numerical  model.  This  comparison 
between  power  generating  experiments  and  numerical  simu¬ 
lations  on  the  same  cell  configuration  has  rarely  been  reported 
before.  The  numerical  distributions  of  current  density  under 
various  operating  conditions  agreed  well  with  the  measured 
distributions  by  segmented  electrodes,  which  had  grooves 
for  the  hydrogen/oxygen  supply  and  were  molded  in  our 
test  cell  being  electrically  insulated.  Hydrogen/oxygen  con¬ 
centration  changes  measured  by  gas  chromatography  along 
the  gas  supply  grooves  also  gave  the  experimental  current 
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Nomenclature 

A 

membrane  active  area  (cm2) 

a 

relative  humidity 

C 

gas  concentration  (g  cm-3) 

Di 

diffusion  coefficient  of  species  i  (m2  s-1) 

d\ 

thickness  of  layer  1  (cm) 

E° 

standard  electromotive  force  (V) 

F 

Faraday  constant  (C  mol-1) 

f 

effective  porosity  of  diffusion  electrode 

Hm-T 

enthalpy  of  species  m  at  T  (K)  (J  mol-1) 

h{  v) 

mass  transfer  coefficient  at  flow  velocity  v 
(cms-1) 

Ijn 

averaged  current  density  (A  cm-2) 

i 

local  current  density  (A  cm-2) 

k 

heat  conductivity  (W  cm-1  K-1) 

MX,M'X 

vapor  flow  rate  (g  s-1) 

Me 

vapor  flow  rate  of  generated  water  (g  s-1) 

Mo$M 

vapor  flow  rate  of  electro-osmosis  (g  s-1) 

N 

molar  velocity  of  v  direction  (mol  cm-2  s-1) 

nd 

electro-osmotic  coefficient 

Pi 

partial  pressure  of  gas  species  i  (atm) 

R 

gas  constant  (J  mol-1  K-1) 

R\ 

gas  diffusive  resistance  of  a  layer  1 

Tmem 

temperature  of  MEA  (K) 

Tr 

transmissivity  of  MEA  (cm2  s-1) 

Fcell 

cell  voltage  (V) 

FNernst 

Nemst  potential  (V) 

V 

supplied  gas  velocity  (ms-1) 

IP  cell 

cell  width  (channel  width  +  rib  width)  (cm) 

Subscripts 

P 

membrane  resistance  (£2  cm2) 

P 

membrane  resistivity  (£2  cm) 

7?  act 

activation  overpotential  (V) 

7?  ohm 

resistive  overpotential  (V) 

Subscripts 

H2 , 02 ,  H2O  hydrogen,  oxygen  and  water,  respectively 

a,  c 

anode  and  cathode 

liq,  vap 

liquid  and  vapor 

mem 

membrane 

MEA 

MEA 

DIF 

diffusion  electrode 

distributions,  which  coincided  almost  with  those  by  the  seg¬ 
mented  electrodes.  Factors  to  correct  the  small  difference 
between  the  measured  and  the  calculated  are  also  discussed 
from  the  standpoint  of  the  physical  meaning  of  the  calculated 
results  considering  factors,  which  are  not  taken  into  account 
in  our  code. 


2.  Measurements  of  membrane  properties  and 
overpotentials 

2.7.  Overview  of  experimental  apparatus 

An  overview  of  the  experimental  apparatus  is  illustrated 
in  Fig.  1  and  is  basically  the  same  as  in  our  previous  paper 
[1].  Fuel  (hydrogen  or  modified  reformed  fuel)  and  oxidant 
(oxygen  or  air)  from  gas  cylinders  were  regulated  by  mass 
flow  controllers,  and  led  to  the  fuel  cell  (FC).  The  humidi¬ 
fiers  placed  before  the  FC,  control  the  dew-point  temperatures 
of  air  and  fuel.  An  Electronic  resistive  load  (PLZ-152WA, 
Kikusui  Elec.  Co.)  was  connected  to  the  FC  instead  of  actual 
loading,  and  measured  the  current  and  voltage  of  the  FC.  The 
MEA  (Gore  Select,  Japan  Gore-Tex  Inc.)  of  20  jim  thick¬ 
ness  is  set  between  the  two  electrodes  and  carbon  cloths 
are  inserted  between  the  MEA  and  separators  as  diffusion 
electrodes.  One  side  of  which  consisted  of  segmented  elec¬ 
trodes  with  grooves  for  air  or  fuel  supply  and  the  segmented 
electrodes  are  molded  in  an  electrically  insulating  plate  for 
measuring  current  density  distributions.  Additionally,  sev¬ 
eral  holes  are  drilled  in  the  grooves  to  measure  concentra¬ 
tion  changes  of  hydrogen  and  oxygen  to  get  another  current 
density  distribution.  More  information  on  the  experimental 
conditions  will  be  given  in  each  section. 

2.2.  Water-vapor  diffusivity  through  diffusion  electrodes, 
transmissivity  and  electro-osmotic  coefficient  through 
the  MEA 

To  obtain  the  diffusive  properties,  the  cell  outlet  average 
humidities  at  both  anode  and  cathode  were  measured  by  sup¬ 
plying  different  humidities  to  both  FC  inlets.  The  amounts 
of  vapor  transport  through  each  layer  Mx  can  be  written  as 
follows,  while  considering  the  mass  transfer  coefficient  h{v) 
between  gas  diffusion  electrode  (GDL)  and  channel  flow. 


=  Aa(u)(Ca  -  Ca,DIFM 


(1) 
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Fig.  1.  Experimental  apparatus. 
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Diffusion  electrode  MEA  Diffusion  electrode 


Fig.  2.  Transport  model  of  water- vapor  in  PEFC. 
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AjIF^a.DIF  ~  Cl.  MEaM 
fifolF 
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Fig.  4.  Change  of  electro-osmotic  coefficient  by  temperature  and  relative 
humidity. 


Mx 


Tr(Ca,  mea  Q,meaM 
^mea 


Mx 


-Pdif(Q,mea  Cc,difM 
Jdif 


(3) 

(4) 


Mx  —  /zc(p)(Cc?  dif  Cc)A 


(5) 


Ca  —  Cc  _  1  1  Jdif  </viea 

Mx  ~  hjv)  •  A  +  hc(v)A  +  DdifA  +  Tr 

=  2/?iO)  +  27?dif  +  ^mea  (6) 

Here  di,  £>dif,  Tr  and  Q  denotes  the  thickness  of  layer  1,  dif¬ 
fusion  coefficient  of  GDL,  transmissivity  through  MEA  and 
concentration  of  vapor  at  boundary  of  two  layers.  The  pat¬ 
tern  diagram  of  this  water  (vapor)  transport  model  is  depicted 
in  Fig.  2.  Considering  the  values  of  Mx  of  Eqs.  (l)-(4)  are 
equal  in  steady  flow,  Eq.  (5)  is  derived.  The  h(v)  increased 
as  gas  velocity  v  increasing,  and  the  diffusive  resistance 
Rh(v)  =  l/h(v)  •  A  was  relatively  much  smaller  than  7?dif 
and  Rmea,  so  that  the  effect  of  the  mass  transfer  can  be 
ignored  in  our  experimental  conditions. 

Measured  water  vapor  diffusivity  for  diffusion  electrode 
Ddif  at  cell  temperature  rmem  =  50  and  60  °C  are  plotted  in 
Fig.  3.  The  diffusivity  did  not  change  against  relative  humid¬ 


Fig.  3.  Change  of  diffusion  coefficient  Ddif  and  transmissivity  Tr  by  tem¬ 
perature  rmem  and  relative  humidity. 


ity.  This  result  is  consistent  with  a  simple  following  equation: 


£>DIF  =  J®air-H20  (7) 

where  £)air-H20  is  water  vapor  diffusivity  in  air,  and  /  the 
effective  porosity  of  diffusion  electrodes.  In  the  case  of  out 
carbon  cloth, /is  calculated  to  be  about  0.25. 

Transmissivity  through  MEA  Tr  is  also  shown  in  Fig.  3. 
Our  measured  Tr  was  about  7  times  higher  than  the  ones 
by  Nguyen  and  White  [2]  and  lower  than  ones  by  Yamada 
and  Morimoto  [3].  However,  the  tendency  of  Tr,  which  rose 
as  relative  humidity  increase,  was  consistent  with  those  by 
Nguyen  et  al.  So  we  approximate  Tr  as  to  7  times  greater  than 
those  of  Nguyen  et  al. 

The  electro-osmotic  coefficient  rid  is  obtained  from  the 
equations  shown  below  and  the  measured  £>dif  and  Tr. 


K 

=  ha(v)(Ca  —  Ca,DIFM 

(8) 

M' 

_  £>DIF(Ca,DIF  _  Ca,MEAM 

(9) 

X 

dp>  IF 

M'x 

Tr(Ca?MEA  —  Cc,mea)^ 

^MEA  +  4^0SM 

(10) 

M' 

_  £>dif(Cc,mea  “  Cc,difM 

(11) 

X 

Jdif  —  Me 

M'x 

=  hc(v)(Cc,Di¥  —  Cc)A  —  Me 

(12) 

Me  and  Mosm  are  water  flow  rate  by  generated  water  and 
electro-osmosis.  Then  the  electro-osmotic  coefficient  rid  is 
expressed  as: 


nd 


4/dsm  F 

18  7 


(13) 


F  is  Faraday  constant  and  I  the  current  through  the  cell. 
Fig.  4  shows  the  rid  as  a  function  of  relative  humidity.  The 
electro-osmotic  coefficient  rid  became  large  when  the  humid¬ 
ity  increases,  rid  at  MEA  temperature  of  50  °C  agreed  with 
that  by  Nguyen  et  al.  rid  by  Yamada  who  performed  similar 
experiment  with  us  was  almost  same  as  ours,  especially  in 
a  case  of  using  carbon  cloth,  however,  the  tendency  of  rid 
against  humidity  was  inverse. 
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Fig.  5.  Change  of  membrane  resistance  by  temperature  and  relative  humid¬ 
ity. 
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2.3.  Ionic  resistivity  and  activation  overpotential  of  the 
MEA 


Resistive  and  active  overpotentials  are  important  factors 
in  PEFC  power  generation  performance.  The  ionic  resis¬ 
tance  was  measured  by  AC  impedance  meter  (SOLARTRON 
SI1280B),  supplying  air  with  the  same  humidity  to  both 
anode  and  cathode  to  avoid  the  movement  of  water  through 
the  MEA.  The  parameters  for  MEA  resistance  are  the  mem¬ 
brane  temperature  Tmem  and  the  flow  relative  humidity  a.  The 
measured  MEA  resistivity  is  plotted  against  relative  humidity 
a  in  Fig.  5.  The  resistivity  did  not  depend  on  the  membrane 
temperature,  but  depended  only  on  the  relative  humidity  a. 
Therefore,  the  following  experimental  ionic  resistivity  p,  a 
function  of  only  relative  humidity  a ,  was  employed  in  the 
following  numerical  analysis. 


p  =  57.3  +  383  exp 


a-  0.089  \ 
0.049  ) 


+244  exp 


/  a  —  0.089  \ 

V  0.268  ) 


(14) 


For  measuring  activation  overpotential  ijact,  the  relatively 
high  flow  rates  of  600  cc  min-1  were  provided  to  hold  the 
low  utilization  rate  of  H2  and  O2  in  order  to  prevent  varia¬ 
tion  of  rj act  along  the  flow  by  the  distributed  current  density. 
?7act  was  derived  by  subtracting  the  cell  voltage  Vceii  and  the 
resistive  overpotential  ir  from  the  Nemst  potential  VNernst  as 
follows: 


had  —  VNernst  ^cell  (15) 

Here,  the  spatially  averaged  partial  pressures  of  active  mate¬ 
rials  were  used  to  determine  VNernst  considering  water  vapor 
pressures.  The  experimental  conditions  for  rj act  measurement 
were  as  follows;  FC  temperature:  60  °C,  cathode  O2  par¬ 
tial  pressure  Po2  \  0.05,  0.2,  1.0;  anode  H2  partial  pressure 
Ph2:  1-0.  Fig.  6  shows  the  change  of  activation  overpotential 
by  oxygen  concentration  and  current  density.  The  activation 
overpotential  increased  almost  linearly  except  for  the  region 
of  very  low  current  density  where  is  not  so  important  for 


Fig.  6.  Change  of  overpotential  by  oxygen  concentration  and  current  density. 

practical  applications.  Therefore,  we  employed  the  follow¬ 
ing  empirical  formula  in  the  following  numerical  analysis. 

0.065 

*7  act  =  -^08T(  +  °'455  (16) 

fo2 


3.  Measurement  of  current  distributions 

Current  distributions  were  measured  by  two  different 
methods  to  confirm  the  reliability  of  the  measurements.  The 
first  method  used  the  segmented  electrodes  cell  to  mea¬ 
sure  the  current  distribution  by  shunt  resistances.  Another 
method  used  gas  chromatography  to  measure  gas  compo¬ 
sition  changes  to  be  converted  to  the  current  distribution 
through  the  decreased  flow  rate  of  H2  or  O2.  N2  mixture 
fuel  and  air  were  employed  in  this  measurement  to  keep  the 
constant  flow  rate  of  N2  as  a  reference.  H2  gas  mixed  with 
20%  N2  (not  CO2)  was  used  as  the  fuel,  because  CO2  is  easy 
to  be  solved  in  water  and  not  giving  the  correct  reference. 
The  cell  configuration  for  current  distribution  measurement 
is  illustrated  in  Fig.  7. 
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Fig.  7.  Configuration  of  cell  and  segmented  electrodes. 


64 


T.  Araki  et  al.  /  Journal  of  Power  Sources  152  (2005)  60-66 


Fig.  8.  Comparison  of  measured  current  distributions  (im=  0.3  A  cm  2). 


5.  The  total  gas  pressure  is  constant,  neglecting  the  pressure 
drop  along  the  flow  channel. 


The  mass  conservation  equations  are  as  follows  by  using 
mole  flow  rate  Mj  of  each  chemical  species  j,  and  current 
density  i  through  MEA.  The  subscripts  a ,  c  and  k  denote  at 
anode,  cathode  and  either  anode  or  cathode,  respectively,  d 
the  height  (thickness)  of  channel  (layer)  and  h  the  channel 
width. 


Hydrogen  : 


dMH2 

dx 


hi 

IF 


Oxygen  : 


dMp2 

dx 


hi 
4 F 


(17) 

(18) 


Water  liquid  : 


d  M\ 


liquid 

water,/: 


dx 


h ^channel 


R(Tk  +  273.15) 


M, 


vapor 

water,/: 


M, 


vapor 

water,/: 


+  Mr2 


P  —  P 
1  ±  \ 


saturated 


water,/: 


(19) 


Water  vapor  (anode) : 


djCga 

dx 


dewater,  a  dMvap0r 


dx 


dx 


(20) 


Fig.  9.  Comparison  of  measured  current  distributions  (im  =0.4  A  cm  2). 


Figs.  8  and  9  shows  the  measured  current  density  dis¬ 
tributions  for  two  average  current  densities  im  of  0.3  and 
0.4  A  cm-2,  respectively.  The  current  density  estimated  by 
gas  chromatography  fluctuated  to  some  extent,  but  was  com¬ 
patible  with  that  measured  by  the  segmented  electrodes  cell. 
Thus  we  could  confirm  the  measurement  reliability  of  current 
density  distribution  by  two  methods. 


4.  Numerical  model  of  PEMFC 

4.1.  Numerical  procedure 

Our  model  basically  consists  of  steady- state  two- 
dimensional  mass,  charge  and  energy  conservation  equations 
with  an  equivalent  electric  circuit  of  the  PEMFC.  The  follow¬ 
ing  assumptions  are  adapted  to  derive  the  governing  equations 
[4]. 

1.  Fuel  cell  is  a  parallel  flow  type. 

2.  Gas  flow  along  a  channel  is  plug  flow. 

3.  Owing  to  the  constant- temperature  water  circulation 
through  separators,  the  separator  temperature  is  constant. 

4.  The  volume  of  the  condensed  water  is  so  small  that  the 
effect  of  condensed  water  at  the  cathode  can  be  neglected. 


Water  vapor  (cathode)  : 

dewater, c  _  _  dewater, C  .  ^vapor  hi 

dx  dx  dx  2  F  1 2 3 4 

Here  Mva por  is  the  water  vapor  mole  flow  rate  from  anode  to 

cathode. 

dA/yapor 

dx 

Rmea  (— Tr(dCu;/dy)  +  itd(i/F ))  -  /?dif,cO'/2E) 

=  - h 

Rmea  +  ^DiF,a  +  ^dif,c 

(22) 


Five  energy  conservation  different  equations  for  five  layers 
(two  gas  channels,  two  diffusion  electrodes  and  one  MEA) 
are  derived,  for  instance,  equations  of  MEA  and  diffusion 
electrode  are  shown  below. 

d  (  drmem  \ 

MEA:  U-  Umem^P  ) 

_  7  Tmem  —  TbiP  a  |  7  Tmem  —  7p)ip  c 

—  ^diff,a-mem  7  r  ^dif,c-mem  7 

^•^DIF,a-mem  ^-^DIF,c-mem 


Nj.aHj-TdiU 


J2n^hj- 


Vh2  A//LHV  "1“  ^Nei  nst  *?ohm  ^aet) 


(23) 
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Fig.  10.  Equivalent  electric  circuits  for  the  PEFC. 
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H. 
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M\ 


liquid 

water,/: 


dx 


+  U(TS  -  Tk) 


(24) 


Here  Nj  and  Hj  represent  the  mole  flux  along  ^-direction  and 
enthalpy  of  the  j  species,  respectively.  U  is  heat  transfer  coef¬ 
ficient.  ?7ohm  is  resistive  overpotential  by  Eq.  (25)  as  a  function 
of  the  ME  A  resistivity  p. 


7?  ohm  — 


(25) 


The  equivalent  electric  circuit  for  PEMFC  is  shown  in  Fig.  10. 
Due  to  the  high  electric  conductivity  of  the  separators,  the  cell 
voltage  along  channel  can  be  kept  constant.  Control  volume 
method  and  tri-diagonal  matrix  algorithm  (TDMA)  were  used 
to  solve  simultaneously  Eqs.  (17)-(24)  with  an  equivalent 
circuit  [5]. 


Fig.  11.  Comparison  between  the  measured  and  the  numerical  current  dis¬ 
tributions  (air  flow  rate:  91  ccmin-1). 
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Fig.  12.  Comparison  between  the  measured  and  the  numerical  current  dis¬ 
tributions  (air  flow  rate:  136  cc  min-1). 


4.2.  Comparison  of  the  current  distributions  with  the 
measured 

The  current  density  distributions  calculated  by  the  numer¬ 
ical  model  are  shown  in  Figs.  11-13.  The  conditions  for 


the  analysis  and  experiments  are  summarized  in  Table  1. 
The  cathode  flow  was  not  humidified  to  observe  the  change 
of  current  density  clearly.  Distributions  of  current  densities 
were  increased  from  the  inlet  as  the  membrane  water  content 


Table  1 


Experimental  conditions  used  in  the  numerical  model 


Figure  number 

11 

12 

13 

Cell  temperature  (°C) 

60 

Anode  flow  rate  (cc  min-1) 

48 

Cathode  flow  rate  (ccmin-1) 

91 

136 

182 

Anode  dew  temperature  (°C) 

60 

Cathode  dew  temperature  (°C) 

Not  humidified 

Fuel  and  oxygen  utilization  ratio  (im  =0.2  A  cm-2) 

0.4,  0.4 

0.4,  0.26 

0.4,  0.2 

Fuel  and  oxygen  utilization  ratio  (im  =  0.3  A  cm-2) 

0.6,  0.6 

0.6,  0.4 

0.6,  0.3 

Fuel  and  oxygen  utilization  ratio  ( im  =  0.4  A  cm-2) 

0.6,  0.6 

0.8,  0.53 

0.8,  0.4 

66 
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Channel  length/cm 

Fig.  13.  Comparison  between  the  measured  and  the  numerical  current  dis¬ 
tributions  (air  flow  rate:  182  cc  min-1). 

increases  with  generated  water.  Then,  they  decreased  from 
around  a  middle  of  a  flow  channel,  since  the  oxygen  par¬ 
tial  pressure  decreased  and  the  activation  overpotential  r/nct 
increased  along  the  channel.  The  position  of  the  current  peak 
was  shifted  downstream  as  the  averaged  current  decreased, 
because  lower  im  means  lower  water  generation  speed  and 
lower  oxygen  consumption  which  leads  to  a  smaller  activa¬ 
tion  overpotential  rjact. 

The  good  agreement  between  the  numerical  and  measured 
results  indicates  that  our  numerical  model  using  the  mea¬ 
sured  water  management  and  power  loss  factors  can  provide 
a  useful  tool  to  estimate  the  PEMFC  power  generation  perfor¬ 
mance.  The  measured  results  showed  a  higher  current  density 
at  the  inlet,  and  a  lower  current  density  at  the  outlet  than  the 
calculated  results,  especially  for  the  higher  oxygen  utilization 
ratio  of  80%.  This  discrepancy  might  be  due  to  flooding  by 
generated  water  at  a  downstream  cathode,  hindering  the  site 
for  electrochemical  reaction.  Further  studies  are  needed  on 
the  empirical  formula  for  the  activation  overpotential,  water 
management  factors  and  diffusion  overpotential  through  the 
diffusion  electrode. 


5.  Conclusions 

The  water  management  factors  such  as  transmissivity  and 
electro-osmotic  coefficient  of  water  vapor  through  the  mem¬ 
brane  electrode  assembly,  and  power  loss  factors  such  as 
activative  and  resistive  overpotentials  have  been  measured. 
These  factors  were  used  to  analyze  our  experimental  results 
for  PEFC  power  generation  tests  by  our  two-dimensional 
simulation  code.  The  code  considers  simultaneously  the 
mass,  charge  and  energy  conservation  equations  with  the 
equivalent  electric-circuit  for  PEFC  to  give  a  numerical 
distribution  of  the  hydrogen/oxygen  concentrations,  cur¬ 
rent  density,  and  gas/cell-component  temperatures  along  gas 
flow  direction.  The  calculated  distributions  of  current  den¬ 
sity  under  various  operating  conditions  have  agreed  well 
with  the  measured  distributions  at  the  segmented-electrodes 
cell.  Hydrogen/oxygen  concentration  changes  measured  by 
gas  chromatography  along  the  gas  flow  direction  have  also 
given  the  experimental  current  distributions,  which  coin¬ 
cide  with  that  from  the  segmented-electrodes.  Degradation 
factors  for  cell  performance  were  also  discussed  from  the 
numerical  results  using  the  simulation  code.  Also  further 
improvements  for  our  experiment  and  analysis  have  been 
discussed. 
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